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A cylindrical solar water heater has been thbricated, 
theoretically analysed and experiaented upon in the present 
study. Hie heater consists of two glass tibes of 15.0 and 
17.8 cm diameter enclosing a G.I. pipe of 11 cm'*, diameter# 
later is made to flow in the G.I. tube whose outer surface is 
painted black, which acts as a solar energy absorber. 

Experiments have been carried cut to test the perfor- 
mance of liie water heater under two different modes of operation, 
namely 

a) natural convection 

b) forced flow condition for mass flow rates of 
50, 90, 150 and 210 kgA» 

Hie day-long collection efficiency under the first mode 
has been found to be 31.5^ vhile in the second mode it lies bet- 
ween 62 - 77^. The maximum efficiency of 76.4^ has been attained 
for mass flow rate of 150 kgA* 

An approximate theoretical analysis has also been carried 
out on the heater for forced flow conditions. The experimental 
results are compared with tiaose obtained ttieoretically, and the 
comparison shows a reasonably good degree of agreement. 
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INTROnJCTION 

1.1 FEED im SOIAE MERGY 
UTIlia^TIDN 

With the increase in industrialization and advancement 
of technology there has been an evergrowing demand for energy 
boih in the developed and the underdeveloped countries, lo get 
some idea of the rate of consumption of energy we shall quantify 
power from 33,000 million tons of soft coal by ’Q' MT . Until 
1850, man had used 9Q. In the next 50 years man had used anotaer 
4Q, The present rate of consumption averages to about 10Q per 
century. By 1‘980 we shall be using IQ every five years. This 
means roughly speaking, that half of all the energy consumed by 
man in the past 2000 years has been consumed in the last cne 
hundred. 

Together with this catastrophe of dwindling fossil 
hiels, is the tempering of ecological balance of nature because 
of pollution udiich comes in the wake of consumption of fossil 
fuels. Thus there is need to look fbr alternative sources of 
Inexhaustible and non polluting energy . 

Ammg the alternative energy resources that are being 
seriousisr; considered at present are t i) Energy from fast 
breeder reactors, ii) Geothermal aaergy and iii) Solar Energy. 

It is believed £2j 1hat if -the fast breeder reactor tedmology 
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comes 'through, then the existing 'uranium resources of 1±ie world 
may last another 1,000 years from now. But this has the serious 
disad'vantage of ecological pollution and scientists have still 
not heaa able to dispose off radioactive waste safely. Oeo- 
thermal Energy also offers an alternative to the Energy cxrisis 
but at best the stored thermal energy in the -world's major 
geothermal areas would, if con^iletely exploited, provide mergy 
for about 50 years. Even In that time the contribution to the 
world's annual consumption would be less than that from tidal 
power which if fUlly exploited is about of the world's 
Energy consumption of 1970. 

Thus we come to the 3rd alternative i.e. the Solar 
Energy. It has its own drawbacks like capital cost of equip- 
ment etc. but from availability point of view it is inexhaustible 
and T/diat is more - non polluting. Thus it becximes quite obvious 
that "this Energy he tapped and used for futur e power needs, 

1.2 AVAIIABIlITy OP SOIAE MERGY 

The amount of Solar Energy incident outside -the atmos- 

i A r —1 

phere amounts -to 1.5 x 10* kwh/year i 3, • But very little 

- 

of this Uaergy is available for use by man. Nearly a third of 
it is lost by reflecticai to space. Another large fraction is 
intereepted by the meteorological thermodynamic engine and its 
associated cl C3Ud belts and cloudy regions. Part of it is 
incident on 70^ of earth' s airface covered by water. Not all 
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of yfaat is left of the solar energy is a'vailahle for buman use. 

1S 

still the aaount of energy falling is about 0.2 x 10 k#i/year 
whida is almost 16,700 times the amount of energy than we now 
consume fi] . 

There are two obvious problems in efficiently harness- 
ing solar energy. Pirst, it is not available continuously and 
thus it necessitates the use of some storage device. Secondly, 
the energy is diffused in form. VThile the amount of energy 
available is enormous, lhat available at any point is not large 
enou^ (solar constant being 2 Cal/aa -min.) to be useful. Thus 
in order that this energy be of use to man it must be efficiently 
and economically harnessed, 

1.3 HAEfflESSIHG Of SOLftR MIEGY 

Primarily Miere are three ways of tapping solar energy, 
depending upon the form in viiich we v®nt the output. They are s 
i) Direct Energy Conversion; ii) Can centra tinn of Solar Energy and 
iii) Collection of Solar Bacsrgy, The first one is used for 
directly getting electric power. Its primary use has been in 
space stations and satellites, ihe major drawback in such devices 
has been their low efficiencies like 10 - 15^ and very high go at 
(goi the order of Rs. 4 3akhs/kw). The second method, namely the 
concentration of solar energy is used for getting hi^ tempera- 
tures e.g, for goaerating steam to xtin pumps or generate power 
or to use it for studying high temperature technolcgy, Presenel 
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lenses and parabolic concentrators are normally used for this 
purpose. The major dra^^ack in Ihese types of concentrators is 
the difficulty in making them, Iheir high cost and the rather 
complicated system needed to make them track the sun. The third 
method, viz, the collection of solar energy is ■the cheapest aid 
■the easiest way of ■tapping this energy, particularly where ■the 
requirement of tempera-tures is of ’the order of 100 ®C, Apart 
from •their usage in i^ace heating or cooling, water distillation 
and fruit drying, -their maximum use has been for solar water 
heating. This can be used in western countries and in some 
northern regions of India in winter. Apart frora this usage, the 
development of -these collectors is useful for a country like 
India from the point of view of making an efficient de-vice for 
generating steam or e-vaporating a liquid of low boiling point 
vdiich in -turn -will drive a pump. Recently Mienel and Mienel 
have found that it is feasible as well as ' economical' to use 
flat plate selectively absorbing collectors -to generate power on 

a large scale. They envisage the to-tal collecting area of II8 3R.; 

2 ' 6 

118 km which will generate 10 megawatts. Secondly -these sort 
of collec-tors can be used as generators for absorption type of 
refrigeration cycle. Thus -to s-fcart wi-th there is a need to deve- 
lop efficiaat solar water heater. 

1.4 REVIEW OP EREVIOOS WCEK 

I Since maximum work has been done in solar water heat*- 
ing, -the liters -ture abounds ih differoi types of heatersy/ These 
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heaters differ from one another on minor points but they can 
broadly be classified into three categories. These are : 

i) Tube type h eater ; 

This is the most common type of heater. Itfe cross 
section is shown in Pig, 1,,1 



Pig. 1.1 1.2 

Yellot and Sobotka f 6l have investigated the perfor- 


mance of this type of solar water heater. The steel pipes are 
simply pressed Into ihenaal contact wiHi the absorbing plate. 

The absorbing plate was coated with a selectively absorbing 
coating as developed by Tabor f?] . There were 2 glass covers. 

The collector area was 2.8 m and steel pipes were of 12 m.m. 
diameter. Baey have reported an efficiency of 38.5^ for the 
ytoole day (natural convection), Por flow rates of 44,5 hg/ii and 
84 hg/li the average daily collection efficiency v^s 44.3f^^iii 
45*5^ respectively. 
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ii) Coriugated sheet type heater 

Ihe heater, Pig. 1.2, consists of a cor lugated sheet, 

coated with black paint, attached to a flat sheet by means of 

rivets.. Water flows in the ’pipes’ foimed by tbe ^eets. Eao 

and Suri sj have euqierimented with this type of heater. Por 

2 

a collector surface area of 2.32 m and the tank capacity of 136 
litres the rise in temperature of water in winter at Roorkee, in 
natural circulation, v«s from 6 °C to 50 (max.) . Ihe efficioacgr 
was found to be The major problem in suc±i type of heaters 

is the leakage from the riveted joints. 

iii) Heater cum storage type 

These types of heaters do away with the storage tanks. 

The heater acts as both. Thus the cost is considerably reduced 

in this case. (Hinappa and Snanalingam ^ 9 J have lasde and tested 

one such heater.; The heater consists of a square coil of 7.55 emV 

diameter and 15.5 m. long in a woodai box with insulation at J 

2 

bottcm and 2 glass covers. Olass surface area is 1.86 m . By 
drawing off v?ater \Bhea its tanperalure reaches 49 °C, it is possi- 
ble to obtain 136 - 230 litres of water/day, Bie collector 
efficienoy has been f>und to be 46^. 

Another type of built-in storage water heater has been 
built and tested at Cmtrai ^id Zone Eesearth Institute, Jodhpur 
by Garg and Rrislmaa j~1ol h* ISie feaa.#ear cra^ists of a galvanised 
iron (20 ^suge) rectan^ijar ’teo^te dljamsions 112 x 80 x 10 cm 
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with a capacity of about 90 litres. This tank is placed in a 
rectan^lar mild sheet tray and insulated from sides and back. 

In front is ordinary glass. The front face of the tank is 
blackoied by lamp black paint. Water temperatures of 55 °C have 
been reported for winter. The efficiency of this type of heater 
has been reported as 75^. 

Chauhan ill j has also built a Imost the same "type of 
heater as above. The capa.city of this heater is 70 litres. The 
maximum temperature attained in the month of April has been 
86 ®0 at about 3.50 p.m. The collector efficiency of 72^ has been 
reported in the forced flow conditions with the mass flow rate of 
76 kg/hr. In the natural convection mode the collector efficiency 
has been found to be between 50 - 58^. 

1.5 MOT IVATIOU FOR THE HIESMT TORE 

There are certain disadvantages with the types of solar 
heaters cited above. One of the important disadvantages is the 
flat surface of the collector itself. Set at a specified angle 
of inclination, the surface receives an effective sunbeam, the 
width of which varies between zero and that of the collector 
during the day, as the axn's position in the sky changes. In 
addition, there exists many avenues for energy losses due to •. 

i) conducticai from the back and edges of plate, 
requiring heavy Insulation for minimization, 


k 


8 


ii) Bafiiaticsa and convection from the flat - plate 
surface and the -water tubes, 

iii) Reflection from and absorption by glass covers, 
caused by -variation in the angle of solar inci- 
dence, and 

iv) ^ionduction -tiirough the insulation surrounding the 
pipes connecting -the collector and the -v/ater tank. 

As a result of these losses, -the average yearly collec- 
tor efficiency is between 30 - 40^. In order to improve -the 
collector efficiency, a cylindrical solar water heater, using a 
cylindrical surface has been fabricated and investigated by 
Vincze 12 | , The collector is a -tube painted black and covered 

by 2 glass tubes with 2,5 cm space between them. The capacity 

2 

of heater is 22.2 1 and the projected area is 0.186 m . Ihen 
placed at an angle of 63 ® to -the horizontal toaperatures as 
high as 58 °C have been reached during the^mnn-th of August. 

Vincze has r^orted efficiencies of the order of 216^, The rea- 
son for such hi^ efficiaacy is given in Chapter IV. The water 
is in Ihe natural circulation. 

As the efficiency appears "tn be very hi^ it has bem 
thou^t appropriate to investigate this type of heater more 
thoroughly, especially since it appears -to be a very good answer 
for airconditioning devices dri-ven by solar energy. 
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The most important advantage of this type of collector 
is that the effective beam width for incident radiation ranains 
equal to the tube diameter over all the sunshine hours of the day. 
Moreover losses die to conduction at the edges are negligible, 
since the area for axial energy conduction is small compared with 
the tube surface and may be well insulated. 

Thus in the present work a solar water heater of cylin- 
drical type has been fabricated. It consists of a G«I. tube 
of 10.15 cm': bore and 11 cnk outer diameter coated with lamp 
black paint. The tube is surrounded by 2 glass tubes of Pyrex 
gle^ss of diameters 15.0 and 17.8 cms. respectively. The capacity 
of the heater is 22.45 litres. The heated section is 1.2 m 
long. 

The heater has been tested for different mass flow 
rates of water steer ting from 0.5 kg/min to 3«5 kg/min. It has 
also been tested for natural circulation. The difference between 
the inlet and outlet temperatures of the water has been noted for 
the viiole day and efficiencies have been calculated. The effi- 
ciaicies have been found to lie betweoi 62 - 77^ for different 
mass flow rates in Ihe range of experimental data. 

An approximate theoretical analysis on ■ftiis heater has 
also been carried out and the results are compared with those 
obtained experimaatally. 
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THBOEEIICAl MALYSIS 

The cylindrical solar water heater under amlysis, 
is as shown in Pig. 2,1, It consists of a metallic tube, vhida 
has been assumed to be coated with a selectively absorbing 
coating | TJ , placed coaxially between two cylindrical shells 
of glass, which act as energy shields. Each inner tube diame- 
ter is assumed to be 43/2 of that of the tube immediately 
surrounding it, so that most of the radiation transmitted by 
the outer lube may be incident csn the inner tube. The axis of 
the tube is maintained at an appropriate angle of tilt in the 
North - South direction, to effect normal solar - radiation 
incidence at noon time, on the design day . Water is forced to 
flow at a constant rate through the metallic tube by some exter- 
nal means. 

It is clear that the total incident intoasity of ra- 
diation on the tube surface, q, varies wiih time of the day. as 
a result, the tanperature of the metallic tube at any section 
is e3q)ected to "vary circumferentially as a function of time. In 
order to reduce the complications that will be intro diced in the 
aralysis of the ten 5 )oral and circumferential variations in solar 
energy input, it has been assumed ihat the total radiation inci- 
dent on the tube surface remains constant at a value equal to the 
average intensity during the sunshine hcairs, for the place under 










)NAL VIEW OF THE APPAR 
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ccaasideration. This assumption, pemits us to treat the systam 
as aie in steady state. Further, the a-vsrage intensity is 
assumed to be due to a ccmbination of the effects of direct 
solar radiation, diffuse sky radiation and "the reflected ground 
and other long wave lenglh radiations. The equations used in 
this calculation are e3diibited later on. 

Since the transmittivity and absorptivity of glass 
vary with the angle of incidence of the radiation, average values 
of these quantities, assumed to remain constant over the whole 
periphery of the glass shields, have to be used in the analysis. 

With the simplifications listed above, it is possible 
to write down the mergy balance equations for the whole ^stem. 

A few other simplifying assumptions used in the analysis are s 

i) The thermal conductivities of the me tall io tube 
and glass surfaces are sufficiaitly high and the 
tube thicdmesses so anell, that the Ihermal resis- 
tances due to the three walls may be neglected; 

ii) Axial conduction through the fluid and the metallic 
tube is negligible, since the product, Be Pr 
(Be = Beynold's number based on the tube diameter and 
the mean velocity of the flow throu^ the tube, 

Pr = Prandtl number for fluid) is large, being 1280 
for lowest mass flow rate of 0.5 kg/fcin. and the well 
thickness is stsallf 
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iii) Concluction losses from the ends of the tubes are 
negligible, since the diameter of 1he outermost 
glass tube is small compared with the solar water 
heater length and tbe eaads are well insulated, 
and 

iv) The fluid properties are invariant in tempera- 
ture because of the small temperature fange 
involved in the heating process. Moreover, the 
viscous dissipation effects have been neglected 
in the formulation of energy balances because 
of the small velocity gradients involved. 

2.1 ENERGY BALANCE EUQA-TI OHS 

Considering the outer ^ass tube, it is seen that the 
tube receives energy due to irradiation from the surroundings 
and due to convective as well as radiative energy exchanges with 
the inner glass tube. It looses energy to the environment 
through convection and radiation. 

-An energy balance on the outer glass tube Iherefccre 

yields s 

^1 ^ 1-2 ^^1 ” ^ 2 ^ ’’’ '^^ 1-2 ®1 ^^1 “ ^2 ^ \ ^ 

= "a ^2 '*2 - ’^ a ) * £-2 ^'^ 2 * ' '' a '*’ 

vfoere and Dg are respectively the diameters of 1he inner and 
outer glass tubes, 2 overall heat transfer coefficient 
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■between the outer end inner glass tubes and P. „ is the "view 
factor for radiation between the inner and outer glass shields. 
The temperatures of the inner and cuter glass shields, respec- 
tively, are T. ^d T„, while h , is the hea-t transfer coefficient 
between the outer cylinder of glass and the ambient air and T_ 
is the ambient temperature. The term ar , is -the Stefan - 
Boltzman constant. The quantities ^'^2 ^2 respec- 

tively the absorptivity and emissivity of the outer glass surface 
for solar radiation and long wave length radiation. 

Pox inner glass tube, similar considerations of energy 
balance yield : 

®1 ^1-2 ^^1 “ ^2^ "^'^1-2 ®1 ^^l"^ “ ^2^^^ 

= -02 (T^- T^) -f- 

+ rr-J) P ^ (T - t/) (2) 

where and respectively the absorptivity and trans- 

mittivity of the inner and outer glass shields, P„ is the metal 
tube diameter, TJ f is the overall heat transfer coefficient 
between the metal tube and the inner glass tube and ^ 

view factor for radiation between the metal tube surface and 

Ihe inner glass tube. 

In writing Eqs, (l) and (2) inter-reflectxons between 
the inner and outer glass tubes have been considered as small and 
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neglected, U^lect of liiese interreflections tends to increase 
the apparent energy losses fron the metal tube surface, and 
hence gives a conservative estimate for solar heater efficiency. 

Since the liquid flowing into the tube is one of rela- 
tively high Prandtl number (4 Pr <C 6)j it bas been assumed 
that the velocity profile is fully developed at the entiy of 
the heated section, permitting use of the thermal entry length 
solution for arbitrary wall energy flux fl?} • the tube 

surface temperature T we can therefore write the equation : 

s 

Tg (x*)- " Tg = [Pg / (2 kp I jf g <x* - 6 ■)•■ Og d1^ 

° (3) 

inhere, 

2 2 

g (x*> = 4 + exp ( - x* ) / ( V A ) (4) 

nn 

and 

= -Cq - (Tg - Vi ~ 

represQits the effective oiergy flux absorbed by tiie metal tube 
surface, assuming that the absorptivity is unity. In Eqs, (3), 
(4) and (5), I is the fluid temperature at the tube inlet and 

'■.'■'S'.' 

X* is a non dimensionaiized axial coordinate measured from the 
tube inlet, given by the equation 

X* = 2 x / (d^ Re Pr) (6) 

On differentiating eq. (3) with respect to x* and 
using the rule for differentiation under the integral sign |J4| , 
we^obtain,; 
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r 


X* 


(2 / D^) a Tg / d X* =J (d / d x*) g - £ )- 4 ^ >- d^ 

O 

+ g (o)- 4 )- (7) 

If we now aibstilute for g (x* - )- sad 4 x* 4 from 

eqs. (4) and (5) into the above and use the integral mean value 
theorou |^ 15 | y we obtain ; 

(2 k^ / Dg) d / d X* = (Ig ( X *) - X *)) 

" cr- [‘ t/ 4 X *} - ■(. x \i } + Cg jz<i “ 

(ig 4 X *> - 4 X *» - ^ ^ ^ 


vdiere, 


and 


, = 5; [«P ( - x» ) - 1] / (V Aj 

n 


( 8 ) 

(9) 

(10) 


C. = 4 + H 1 / (>* ^ A ) 
n 

— * 

while X , represents a point in the Interval 
0 ■C. X k_ X , at which the various temperatures are to be eva- 
luated. further, in eqs. (4), (9) and (lO), the quantities 
are eigen values and the quantities A^ are constants, listed in 

Li 53 • ■ 

The mixed mean ten^jerature, T^, of the fluid flowing 
throughthetubeisgivenbyequation;|^13f- 

d T^ / d X* i= (2 Dg / kp q^ 4 X* >• (n) 

kj being the fluid iheri!®.! conductivity. 
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If eqs, (2)5 (3), (s) and (11) are now solved simul- 
taneously to determine the temperatures, I and T , as 

functions of x, we can eva-luate the heater efficiency (betv/een 
the inlet and the point under ccaasideration) by writing s 

7 ] = m Cp (T^ix^ - / (b^x q) (l2) 

# 

where m, is the mass flow rate of fluid flowing through the 
tube, and is the entry temperature of the fluid. 

2.2 EViaDUATIQU OF COISIANTS AND HEAT ELUZ q 

Before the solution of the above equations be attempted 
it is necessary to find tlie constants like '^1-2’ ^ i^^l-2 

etc and ttie solar radiation input. We shall calculate each of 
these one by one. i 

Solar Energy Radiation Input q 

The total radiation falling on a horizontal plane at : 

I.I.T. Kanpur has been recorded for evemy day by the Civil i 

Engineering Departmait on the pyraaometer supplied by Indian 
Meteorological Department. Graphs have been taken from them both 
for direct and diffuse radiation. 

The total radiation on the collector consists of; 

i) direct radiation on the top surface, ii) diffuse ^y radia- 

• ^ I: 

tion on the top surface and iii) ground reflected radiation on 

the bottom surface, s 

■ i 

Tirr^ KAMPU^ 

central UBRAlj 

; Acc. No. A 3&6*l 
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The heater is inclined to the horizontal at an angle 

r* 1 

calculated from the equations : j 

Sin f = Cos I Cos S' Cos H + Sin L SinS 

Sin (z/" = Cos S Sin H / Cos , and (l5) 

Cos 0 = Cos p Cos / Sin + Sin |S Cos ^ 

vtoere 1 = latitude arglej for Kanpur it is 26°26', £= declina- ' 
tion angle and H = hotaap angle = 54’ . The sunshine hours are 

taken to be from 7 AM to 5 PM. It has been assumed that the de- 

clination angle changes very little between 29th May to June 12 
(test dates), ranging between 21° and 23°. Thus the average f-ngle 
has been taken to be 22°. 

Assuming that the direct solar radiation should be normal 
to the tube axis at 12 noon for maximum utilization on the design 
day (i.e. 0 = 0 or Cos © = l), the equations above yield; 

’ic 6° 

which is the inclination of the tube to the horizontal. 

Thus the direct normal radiation, as seen from the 
graph of insolation, be miltiplied by l/CosJf' to make it normal 
on the inclined surface. 

The diffuse radiation falling on the inclined surface 
shouldbemultiplied by a factor : j 1?] 

. R = (l + Cos If ) / 2 

S',' ^ 

so as to make it normal on the inclined surface. 
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The ground reflected radiation, falling only on the 
under side of ihe heater, is found as follows : 


R, 


RE 


( E 


m 




E 


JL 

D g 


■where E^ = Direct normal radiation, R^ = normal dif Rise ra- 
diation, = shape factor between the ground and the heater = 

•81 and r = ground reflectance = .12 for gicavel roof 

6 

top, lie] 

The total radiation, for the day of experiment, falling 
on the heater can "thus be calculated as : 

i) Direct radiation s This radiation is assumed to 
fall only on the projected area formed by an arc 
. of 120° of the cylinder. The rest gets all re- 
flected from the surface. This area for the heater 
2 

is O'. 186 m . 

ii) Diffuse radiation is assumed to fall only on ihe 

2 

upper portion and the projected area is 0,214- m 

iii) Oround reflected reflection is assumed to fall only 

on the under side of the collector and liie projec- 

2 

ted area is 0.214 m . 

Overall Heat Transfer Coefficients 

The overall heat transfer ooeffieients, „ and U 

1-2 S—1 

appearing in eqs, (2), (3) and (s), have to be obtained from the 
curves on overall film conductances for free convection flow 
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between ccn centric cylinders provided by Jakob |_19} . Under 
the cmditions of the presait computations where Dg = (2 / 's/3 ) 

= ( 4 / 3 ) D > the overall film coefficients are very nearly equal 

o 

to the respective thermal resistances caused by the presence of 
annular air gaps, we can therefore write for U „ and U ^ the 

/ S'**! 

approximate equations ; 

“1 "l-2 = 2 \ / Hi ^2 / H,). ‘-na = 2 / to (E^ / 

(15) 

Outside Heat Transfer Coefficient 


Assuming the average wind velocity of 12,15 km/h in the 
East-West direction for Kanpur | 20 1 the heat transfer coeffi- 
cioit hg^, between the outside air and outer glass ^ield of dia- 
meter, ©2 can be calculated by using the equation \_2ll 


a 2 


= a ( He )* 


vhere k„ is the thermal conductivity of "the ambient air 

'Q. ' 

and the Reynold's number, Re is based on the outer tube diameter 


and the design wind velocity. The quantities, a and b, are cons- 


tants which depend upon Reynold's number. The value of h^^, cal- 

2 

culated in this manner comes out to be 15 kcal/hr - m - °0. 


Shape Eactors 

The calculations involve shape lectors for a system of 
two concentric and coaxial cylinders, lidiere the inner surface of 
the outer cyliniter radiates to itself. Shape factors for thi^^ 
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case has been presented in the fora of graphs by Hottel and 
Sainfin"- [_2^ and have been used In the present calculations. 

A knowledge of these shape factors enables the deterainastion 
of ^s-l’ ^ Eqs. (2), (5) Sind (8), by ihe 

use of standard eapressions for radiant energy exchange between 
coaxial cylinders with gray inner surfaces f 22j . The values 

of ^s-1 °^^oulated in this way are : 

^1_2 “ 0.602', and = 0.895.- 

The values of transaittivities, absorptivities and 
emissivities of glass and oetal surface (with a selectively 
absorbing blade coating), are listed below in Table 2-1. 

TABLE 2-1 

Transaittivities, Ab scarp tivi ties & Emissivities 



Glass 


Single 

0.88 

0.05 0.06 

0.94 

Double 

0.81 


■iHm' 


Metal Tube 


0.0 1.0 0.91 0,09 


2.3 METHGL OF SOEITIOH 

The coupled and non linear differential equations in- 
volved in the determination of T^, Tg ®^d T^, preclude the 
possibility of obtaining analytical solutions. The simplest mode 
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Of attack is therefore to otrtaia nmerical solutions for para- 
metric -values of the tube diameter, D , -the mass flow rate, m 

s 

and -ihe tube length. Using the da.ta as specified abo-ve along 
with the specified diameters #iich are 11 cms, 15 cm and 17»8 cm 
for metallic tube, inner glass shield and outer glass ^ield 
respectively, and the ekss flow rate of the liquid through the 
tube, an iterative, step-by-step calculation may be carried out. 
The steps in-volved in the computation are the follovTing ; 

a) C!hoosin.g a sufficiently small x, say i cm., and assuming the 
metal surface at the inlet, equal to that of the entering 
water, the simultaneous equrtions for T^ and T^ at the inlet 
are determined numerically. The technique used in -the cal- 
culations is Newton's method of tangents, assuming that the 
initial guesses are sufficiently close to the roots [ 23^ • 

b) Using Bqs.(2) -(8) as well as T^ and T^ as determined above, 

a first estimate of T 4 2^^ is ob-bained. In this calcula- 

s . 

tion it has been assumed -Ihat x is approximately equal to 
x*/2, since x* is small (13, 9 x 10 ^ / ra.) snd there is not 
likely -to be much of difference between x and x*" / 2. 

c) After settir^ T^ 4 / 2 •)■ = T (o) + T_ 4 ^ / 2, 

steps (a) and (b), are repeated until the newly computed 

values of T ^ x* ^ agrees with -that computed early closely. 

S'.', 

The values of T^, and Tg corresponding to T^ 4 x* / 2 and T^ 
^ xXr ^ as well as the liquid tenperature, T^ 4 me* ^ are de- 
termined by using 14a e appropriate equations. 
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d) The efficiency of ihe solar vvater heater in the interval 
0 to X is 

^ °P ( \ 4 3£* } - Te ) / (lig q x) * (l6) 

e) Now, the origin may be translated to the point x = x*/2, so 
that the new starting values for the computetions are 

Tg (O) = Tg 4 xV2 ) , T^ (O) = T^ ( x*/2 ) and T^ (o) = 

^2 ^ ). The calculations at steps (b), (o) and (d) are 

repeated, until the temperatures along the tube surface as 
well as the mixed mean fluid teirperatures, and •&€ efficien- 
cies at all values of x, are obtained. 

The computer programme to eva.luate the temperatures 
using ihe above scheme is listed in Appendix I. 

The computations may be stopped, once the end of the 
tube is reached which is 1,20 m in 1he present set up. However, 
if the purpose of the calculation is to determine the solar heater 
length needed for maximum aaeigy collection when the tube diameter 
eaad the mass flow rate ai-e specified, it is necessary to continue 
■the calculations imtil -the newly computed fluid tempera-tures 

^ differ appreciably from -that one step computed 

earlier. Por •the "tube diameter of 18 cms and mass flow rates of 
0,5, 1,5, 2.5, 3.5s 4.5, 5.5 and 6.5 kg/min, ‘ the optimum length 
has been found [ 24 J . Por mass flow rate of 0.5 kg/min. ‘ the 

optimum length comes out -to be 6 m, the difference between T^ at 
6 m and that at 7 m being less -than 2^, 


* Please see page 24 



In -tile present viork the conputations have been done 
for mass flow rates of 0.5, 1.5-,,c2.5 and 3.5 kg/mn. The 
comparison of results between those computed and found experi- 
mentally are diown in Pig. 4,6, . 


* Hitt «Kpnii»l«a i« m% t«r mtwml mm* fim •INfi«t«wi3 
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afaaaiirw iv 'Va fptrapia'<iii ^pmaa^itiiiw wia . mwm laiaiaR! iBpaaiaapipai> 

f^m fmaimmtimm at aataa at tta mA af Aiir 

1 t i a p i>f*t «»»* af aatar la <>a iiaAag 

f « fatal nriAatliM (iliaat» ilffiMa ani pwMrii laflaatai) 
iaalMI m Hmd taaNir <iaflBg lia iigr* - ' 




CHAPTER III 


EXPERIMEINTAL SET UP 

In order to check the theoiy and also to study cylindrical 
solar water heaters more thoroughly, a heater has been fabricated 
and ejqjerimgnted upon. The line diagram, of the set up is shown an 
5'ig* The main parts of the set up will be described under 

two headings, i) Appaietus and ii) Ins tram aa tat ion. 

3.1 APPAHATUS 

i) G.I. Tube 

A tube 2,5 m long, 10,15 cm. inner diameter, and outer 
diameter 11 cm,, made of Galvanised Iron is used for heating the 
T^ater, In the theoretical analysis of the heater it has been 
assum«i that the flow of water in the tube is fully developed. 

Hence a sufficiently long entry length must be provided at the 
tube inlet to oisure a fully developed flow. 

The hydro dynamical entry laogth can be calculated by 
using the equation s [^25 j 

X = 0,03 Re D 

vJhere x = entiy length, Re = Reynolds number and D = diameter 
of the tube. The mass flow rates for whddh the heater has been 
ejiperimented upon vary from 0,5 kg/min to 3.5 kg/min.. Eor 3.5 
kg/min. mass flow rate, the entry laiglh comes to 3.2 m. Since 
such a long tube was not available the longest tube available, 

Mhich is 2.5 m Icng has been used. 






TO’ ' 




till 

pi»B 
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The heated section of the tube is 1.2 m. Thus the 
ranaining 1.3 m length of the tube acts as ihe entiy length. Ihr 
mass flow rates of 0.5 kg/min. and 1.5 kg/min. the entry length 
as found by equation above comes to 0.45 m and 1.35 m respectively. 
Thus fully developed flow as assumed in the theory can exist only 
for flowrates of 0.5 kg/min. and 1.5 kg/min..,,. 

Since this lengiii is meant only for developing the flow, 
it has been insulated by 5 on thick glass wool so that sun does 
not heat up this portion. By making simple calculations it has 
aiso been seen that the amount of energy input in this section, 
becaise of direct, diffuse and ground reflected radiation, comes 
out to be less than of that incident on the heated section. 

The portion of the tube to be heated has been painted 
black with a paint vhidi is a mixture of lamp black and black 
aiamel paint. 

Water is allowed to flow from Ihe lower end of the tube 
■?4iich is inclined at an angle of 6° to the horizontal, as ^ovn 
in Pig. 3»1. The flow rate has been measured wilh ihe fcalp . f 
of a flow meter. 

ii) Wooden Panels 

These panels are mounted on 1he tube and clamped to it 
by means of 6 bolts, 3 on each one of ihem, as shown in Pig. 2*1. 
At the bade of ihese panels are the thermocouple wire connectors. 
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In order -that iiie connectors are not exposed to sun they have 
been shaded by vdiite packaging material, 

iii) Glass Shields 

The two glass fields are of-a^skD-ian'S. cm-adiamBter 
respectively. They are made of pyrex glass with average trans*-' 
mittivity of 0,89. The glass shields are 2 ram thick. The air 
gap between the G.I. tube and twD glass shields is 2 cm and 1,5 cm 
respectively. These shields are mounted on steps made on the 
wooden panels, as shown in Pig. 2,1. Eubber gaskets have been 
placed at the rim of the panel so that contact between glass 
shield and panel is good, thus preventing any enterance of dust 
or dew. 

5.2 maCEUMENTATIOI 

i) Taaperature Measurements 

The temperatures of vieter as well as the G.I. tube and 
glass Shields have been measured by calibrated 24 gauge insulated 
copper constanton thermocouples, 16 tliermocouples have been 



fig 3.2 
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are fixed on outer and inner glass fields by means of an adhesive} 
epoxy resin. 

The water inlet and outlet temperatures have been mea- 
sured by dipping the liiermocouples in inlet and the outlet pipes 
as shown in Pig, 3-1. 

The ambient temperature has been measured by putting a 
thermometer, with a least count of 0.1 °0, In shade. 

ii) Recorder 

A Leeds and Forthrup millivolt recorder with a least 
count of 0,05 mV has been used to record tie temperatures, 

iii) Solar Radiation Measurements 

d 

Solar radiation has been measured by Epp^ pyranometer 
installed by National Meteorolcgical Laboratory Poona at the 
Civil Bagineering Lepartment of I.I.T, Karpur. The charts frtom 
this recaorder for both direct and diffuse radiation have been 
used to calculate tie amount of radiation falling on the solar 
water heater. The radiation as read from these charts is that 
incident on a horizontal surface. The radiation incident on an 
inclined surface corresponding to the inclination has been cal- 
culated therefrom by using appropriate factors as shown in 
Chapter II. 

3.3 MOLES OP TEST 

The performance of the heater has been tested under two 


cmditions s 
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i) Natural circulation 
ii) teced circulation 
Natural Circulation 

In this mode, the pipe is filled with water in the morn- 
ing and the tap at the lower end of the tube is closed. The 
teiiipcjra tires of the tube and the glass shields are recorded. At 
the end of the day, the inlet and the outlet of the tube are 
connected by a small pump (o.025 hp) through polythene pipes so 
that the water is thoroughly mixed. The final temperature of the 
water is thaa recorded. 

Ihrced Circulation 

The water is forced in the tube from below. There is 
no use of pump as the tap pressure is enou^ to force the water 
to the maximum flow rate of 3.5 litres/min. Ten^sera tuxes both 
for inlet and outlet water as well as those of the tube and the 
glass shields are measured. At the end of the day the water in 
the tube is mixed by the pump in order to measure the total ther- 
mal energy gained by the system. The heater has been tested for 
four mass flow rates, 0.5> 1 . 5 ^ 2,5 and 3,5 kg. /min.. 



GEmim IV 


RBSIJBDS -ill) mSCOSSIGHS 

4,1 HESUBrs 

The experimental phase of this study csonsisted of tests 
conducted during the months of May and June, 1974. luring the 
test period, a wide tange of olioatic conditions were encountered; 
particularly the onset of early monsoon made it difficult to carry 
out tests on all days. Thus only the data chtained on the days 
'fhm there was uninterrupted sunshine is presented in this work. 

The heater 'was tested in two inodes viz. natural convec- 
tion and forced oonveotion. Table 4.1 shows the experimental 
data and results for natural convection. 

Table 4.1 : Experimental Data and Hesults for Test in 
Satural Convection (May 18, 1974). 


Water temperate e] 

°a , ! 

! Time | 

^ ^^1 Inso lat ion 

onglel 

(koal) 

[gained Icximcy 

1 


IFpcxa 1 To 1 

{A.M. 1 P.M.K 

jDireot {Diffuse 5&rcnand {'by % 

? f freflectedl (kcal) 1 

3i.5 

55.3 

7,20 6,10 

8° 1040 487.15 

172.5 535 31.5 


Figure 4.1 shows the plot of the experimental data for 
natural convection. The direct and diffuse solar radiation fall- 
ing on the tilted collector are plotted, together with 1he plate 
temperature and the ambient temperature* The maximum plate tempe- 
mture of 63.5 °C was obtained at 3.00 p.m. after which it started 




Insolation^Kcal /m 



lUXJ 


Time of day, hours. 

' Direct- normal radiation ori collector surface 
'’'^^-Piffuse sky,.. '• ' 

•Plate temper ature,.,';,,,^..x 

4 .1 Ex per f m ental^^ jStts 
circuiation:"%^ '' 


:i' ;Ki3n'(>ur 


y ,18 1974 for natur 


900 


a30 
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to decrease since the heat losses from the plate to Idle surzound- 
iags were more than the input by solar radiation. 

The efficiency obtained in this case is rather small 
being 31*5^. A possible reason is that since the angle of tilt 
of the heater to the horizontal is small (8°)j convection 
currents do not develop readily. Also since the upper surface 
of tile tube is heated there is no radial convection. More vigo- 
rous convection could have been obtained if the angle of tilt 
had been around 30°, but then this would have been at the ex- 
pense of less solar input to the heater, becanse then the solar 
input would not have been normal to the surface of the heater. 

Since the efficiency in natural convection is small 
a better indication of the effectiveness of the heater would be 
to test it under forced convection. This does not mean that 
the efficiency of such a type of heater will be small in natural 
convection, but it simply means that under the cenditions of 
present esperiment viien the inclination eingle is 8°, it will be 
small. Table 4.2 tabulates the results for test in continuous 
flow for different mass flow rates, 

Figures 4,2 - 4.5 ^ow the plot of eoqperimental results 
for forced convection. Together with liie total radiation falling 
on the heater, the outer glass tanperature and the difference 
between the water temperatures at ihe outlet and the inlet are 


plofted 



12,6.74 e.eo 5,40 6° 210 961.8 0,0 965 441. C 151.0 1557 961.8 62.0 
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TABLE 4.2 : Experimental Data and Results far test in Continuous Flow 




o Direct insolation on the surface of collector cp =6^^ 

' •Diffuse insolation » , cp ^s® 

' ;Xj Outer glass temperature - _ ■ - ^ 

:':,d"Ambient temperature ^ 

A Di f fer ence between ^ water, . at outlet & inlet 

r. > 

FIG. 4.2 Fxperimentah're^ults^#;' 5 pf||S^^^^^^^^^^ '29 M'ay,1974 for 

^ ; 'mass, Jlow;rigTe''of 


,uJ-. .v'>^ 




ilii 


Temp. difference f 





e of day, hours 
)eratur <2 

diation on collector surface 4^ = 6^^ 
ntion 5> . ” 4^ 6° 

iture 

en water temperature af outlet & inlet 






rT' 'I , 
. ‘* 5 ) 


Temperature,' 




nsoiation ,Kca 


800 



Time of day, hours 

o Direct .ra.diation on collector 
• Diffuse sky ' ' 

7 Outer glass temperature ; : 

□ : Ambient temperature , 

h ' p'i f f er'e nee , fcefw een. dutlet & dnfet 











o Direct norma! Yadfatton on the collector surface 
• Sky diffuse radiation on the collector surface 
X Outer glass temperature 
D Ambient temperature , . 

A Difference between water temperature dt, outlet 8« Inlet 



df „liqn'pur,jyne 12^1974' forYma; 


Temperature °C 
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In all these figures it is seen that the difference 
between 1he outlet and the inlet temperatures of water at the 
initial stage of the experiment is negative and that the tine in 
which this difference remains negative decreases as the mass 
flow rate increases. A possible explaination for this can be 
given as follows. In the morning the heater is filled with 
water at a particular temperatare. The inlet water is fed in 
from a tap which is connected by a long system of Gr.I, pipes 
to an overhead tarlc. Thus 1iie inlet waiter is warmer than that 
in the heater pipe. In the early sts.ges of the eaperiment, 
this hot incoming water gets mixed v/ith the colder water inside 
the heater and this becomes cooler. Also in the momiiog ihe 
amount of sun’s radiation isnot enough to provide as much energy 
as it has lost in the insulated portion. For the lowest mass 
flo?/ ra.te of 0.5 kg/min. the time during which the outlet and 
inlet water temperature difference remains negative is about 
3 hour. As liie mass flow rate increases the difference between 
temperatures of the incoming water and that in insulated section 
goes on decreasing. Also the time of contact of the incoming 
water w'ith that in the insulated section reduces. Therefore 1316 
time during vshich the difference between the outlet end inlet 
water temperature remains negative, decreases. For mass flow 
rate of 3.5 Isg/min. this time reduces to less than an hour. 

From Table 4.2 it is also seen tinat 1he heater effi- 
ciency increases with the increase in mass flow rate. As Hie 
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mass flow rate is increased the heat transfer coefficient also increases 
and thus there is better transfer of heat from the tube to the fluid* Also 
at higher mass flow rates the temperatujres of the plate are lower than those 
at lower mass flow rates and thus there are smaller losses from the plfttc 
to outside. A combination of these two effects loads to the observed rise 
in efficiency. If a second order regression curve is drawn through the data 
points of efficiency vs. mass flow rate then all the points lie within + 5^ 
of this CTirve. At lower mass flow rates of 0.5 and 1.5 kg/min the trend 
shows an incFoase in efficiency with increa; e in mass flow rates because of 
reasons stated above. For higher mass flow rates the trend is difficult to 
predict since besides the scattering of points due to exiierimental error 
there is the problem ©f flow being not fully developed because of inguffi'" 
cient entry length. 

Figure 4.4 shows the experimental data for a mass flow rate of 
150 kg/h. There is a sudden dip in the incident radiation of the sun at 
aboL^t 1 p.m. This was because of clouds. This effects the water temperature 
difference which also dips at 1,30 p.m. 

Finally the experimental results are compared v^ith those obtained 
theoretically. Figure 4,6 shows the comparison for efficiency, plate 
temperature and difference between water temperature at the outlet and 
inlet, for varying mass flow rates of 0.5 kg/min, to 3.5 kg/min. 

Table 4.3 lists the various quantities, bothi experimental 
and theoretical. The theory compares reasonably well with tbe 
experimental resdlts obtained. The difference between the theore-" 

tical and experimental results for mass flow rates of 30 kg/h to 
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Plate temperature 


Theoretical 


0 Experimental 
Difference in water temp, at outlet & inlet 
— Theoretical v Experimental ,, , 
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150 kg/h lies within 11^. She iiieoretical efficiency is higher 
since we have assumed in 1he analysis 1hat the absorber plate 
is coated ?dth a selectively absorbing paint vhich has absorp- 
tivity of unity and very low emissivity. In the actual test 
the emissivity of the lanp black coating is also nearly equal 
to unity. Hence one may expect higher losses due to radiation 
end a reduced efficiency. Also we have assumed in theoretical 
calculations 1hat the tube temperature is uniform at any given 
section vhei’eas this ca-nnot possibly occur in practise owing 
to non-uniform irradiation over the tube periphery, and the 
relativelj'' small ihemal conductivity of the tube. 


Table 4.3 ! Comparison of Experimental and Theoretica.l 
Results of Heater under ibrced Convection 
Mode. 


Mass flow \ Efficiency io 

rate kg/h | 

feheore'ticall Exptl 


30 

74.5 

67.0 

54.0 

49.5 

4.12 

3.4 

90 

78.4 

70.3 

46.5 

44.7 

1.58 

1.29 

150 

78.35 

76.4 

45.3 

44.0 

1.05 

0.94 

210 

78.3 

62.0 

44.5 

43.5 

0.64 

0.51 


Jpiate Tecp. °C ][l)iff. of outlet and " 

I |inlet water temp. °G/h 

iTheoreticaliExptli Theoretioali l^tl. 


It is worthvhile to comment here on the efficiencies 
reported by Vincze j 12 | . He has reported efficiencies as hi^ 

as 216^. There are two probable reasons for such absurdly hi^ 
efficiencies. One is that Vincze has assumed the radiation falling 
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The theory compares reasonably well with the e3q:.eri- 
mental results. The scf^llest difference of 2.8/S is found for 
mass flov? rate of 210 kg/h vdaereas for mass flow rates of 30 
and 90 kg/la the differences lie in the range of 10 - 11^, 

Still from purely heat transfer point of view there is ample 
scope of improving the theory, and it may be taken up as the 
next stage of investigation. The variation of solar radiation 
with time, together vd-th the circumferential temperature va.ria- 
tion may be incorporated in the refined theory. 

The present study is also a first step tewards making 
an efficiikit solar collector device viiich can be used as a gene- 
rator of a vapor absorption refrigeration system. The tempera- 
tures obtained in the present study are not high and in order 
to obtain higher temperatures a parabolic concentrator can be 
put belov; the cylindrical heater. The cylindrical heater has 
one more advantage over the flat - plate collectors, if used 
for purposes of generator, and that is that the cylinder can 
withstand higher pressures easily while ii flat plate collectors, 
bulging of the absorber plate and consequent leakage may set in, 
if high pressures are used. 
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w 31 AUGUST 7 

IBFTC MAIN 

dimension 

COMMON X,TS,TGI,TGJ,T 
PFAL NETAUZOC) 

READ, C VM( I ) ,AM( I ) ,1=1 
RS=0.5 

31 ■X=(13.9E-5)/RS 

READ,TS,TGItTGJ 


SHGjTGI HC, A F,A0,TS 


AQ = A. 

DO 20 1=1,5 

AP=AP+(EXP< (“ VM( I) )*¥ )•”;.)/ (VM( n* AM ( n ) 
AQ=AQ+1,/(VM( U*AM{n) 

TE = TS 
TM=TS 

CALL ALGEQfTS ,TGI,TGJ) 

TSN = TS 

TGIN=TGI 

TGJN=TGJ 

TSI=TS+G(X,TS,TGI»#rAP+AQ) 

CALL ALGEQ{TSI,TGI,TGJ) 

TSH=0,5*{TSI-TSN|+TSM 

TGIH=0.5*(TGI-TGIN)+TGIN 

TGJH=0.5«(TGJ-TGJNI+TG JN 

TGIH0=TGIH 

TSHO=TSH 

TGJH0=7GJH 

CALL DIFEQ 

TSH=0,5*(TSAN-TSN J+TSN 
CALL ALGEQ<TSH,TGI,TGJI 
TGIH0=TGI 

i,TGJHO«TGJ ‘ ■ 

TS,MD=TSH ^ ' 

I'Jw'l, A,'', y. . 

CAll ,,ALGfiQ { TS H ,TGI ,TG J > ' ' 

TGIH<=TGI- , 

.T64H®TGJ 
G0 TO t3 
TGI'»TGIH0 

T{»4?»TGJM0 ' , ' 


ALGEQCTSAN,TGI»TGJ> 


■GJHC=TGJ 









